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摘要 :【 目 的] 建立 基于 小 波 神经 网 络 病虫害 预测 预报 模型 ,对 提前 采取 防 病 防 忠 措施 减少 农作物 病虫害 损失 、 提 
高 农作物 产量 与 质量 具有 重要 意义 。[【 方 法 ] 本 研究 以 山西 省 运城 市 商城 县 1980 -2014 年 麦 蚜 发 生 程度 和 气象 因 
子 数据 为 基础 ,采用 主 成 分 分 析 法 从 40 个 基础 气象 因子 中 整合 形成 9 个 新 的 自 变量 输入 模型 ,采用 试 凑 法 筛选 隐 
含 层 节点 数 ,用 1980 - 2009 年 的 数据 进行 网 络 训练 ,对 2010 - 2014 年 麦 蚜 发 生 程度 进行 回 测 ,建立 了 以 Morlet 小 
波 函数 为 传递 函数 的 小 波 神经 网 络 模 型 ,并 与 以 Sigmoid 函数 为 传递 函数 的 BP 神经 网 络 模型 进行 了 比较 。【 结 
果 】 小 波 和 BP 神经 网 络 两 种 模型 对 训练 样本 的 平均 拟 合 精度 均 有 10 年 以 上 超过 80% ,两 者 MAPE 值 分 别 为 
89. 83% 和 83.07% ,MSE 值 分 别 为 0.0578 和 0.6192。[【 结论】 两 个 模型 都 能 较 好 地 描述 麦 蚜 发 生 程 度 ;从 预测 精 
度 和 模型 的 稳定 性 来 看 ,小波 神经 网 络 好 于 BP 神经 网 络 。 

关键 词 : 麦 蚜 ; 小 波 神经 网 络 ; BP 神经 网 络 ; 发 生 程度 ; 预测 

中 图 分 类 号 : Q968 ”文献 标识 码 : A ”文章 编号 : 0454-6296(2015 )08-0893-11 

Forecasting model for the occurrence degree of wheat aphids based on 

































































wavelet neural network 

JIN Ran, LI Sheng-Cai ( College of Agriculture, Shanxi Agricultural University, Taigu, Shanxi 030801 ， 
China) 

Abstract: [Aim) This study aims to build up a pest and disease forecast model based on wavelet neural 
network (WNN), so as to provide a basis for taking measures to prevent pests and diseases, reducing 
crop damage by pests and diseases and improving quantity and quality of crop yields. 【 Methods)] Based 
on the occurrence degree of wheat aphids from 1980 to 2014 and the meteorological factors in Ruicheng 
County, Yuncheng City, Shanxi Province, we integrated and created 9 new independent variable input 
models from 40 fundamental meteorological factors through Principal Component Analysis (PCA) and 
screened hidden layer nodes by trial and error method, conducted training with data from 1980 to 2009 
and retested the occurrence degree of wheat aphids from 2010 to 2014. Finally, the study built up a 
WNN model by taking wavelet function as transfer function and contrasted itself with BP neural network 
( BPNN) model which takes Sigmoid function as transfer function. 【Results】The average fitting accuracy 
of both models, namely, WNN and BPNN models, were above 80% in at least 10 years. Their MAPE 
values were 89. 83% and 83.07%, and their MSE values were 0. 0578 and 0. 6192, respectively. 
【 Conclusion 】Both models can well illustrate the occurrence degree of wheat aphids. As for the forecast 
accuracy and model stability, however, WNN is better than BPNN. 
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多 年 来 ,昆虫 学 家 和 植保 工作 者 致力 于 病虫害 ”统计 预测 法 等 传统 预测 预报 方法 ( 张 孝 闵 等 ， 
测报 方法 的 研究 ,发 展 了 经 验 预 测 法 、 实 验 预测 法 和 1985) 。 但 由 于 害虫 所 在 的 生态 系统 具有 复杂 的 非 
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线性 特点 , 易 受 环境 因子 及 害虫 自身 生长 发 育 和 天 
政 发 生 情况 等 多 种 因素 影响 ,传统 的 线性 预测 预报 
方法 很 难 达到 理想 的 效果 。 于 是 ,一 些 专 家 ,学 者 开 
始 将 现代 非 线性 理论 运用 到 害虫 测报 领域 ,将 传统 
的 动力 学 理论 .数理 统计 与 现代 计算 技术 相 结合 ,发 
展 了 人 工 神经 网 络 、 相 空间 重 构 预测 法 、 小 波 分 析 、 
支持 向 量 机 等 病虫害 测报 新 方法 ( 马 飞 等 , 2002a， 
2002b; 袁 哲 明 等 , 2008; 张 永生 , 2009 ) 。 

人 工 神 经 网 络 (artificial neural network，ANN ) 
模拟 人 脑 结构 设计 ,是 人 脑 的 一 种 物理 抽象 .简化 和 
模拟 。 它 采用 并 行 分 布 式 系统 ,具有 良好 的 自学 习 、 
自 组 织 、 自 适应 特性 以 及 处 理 非 线 性 问题 的 能 力 , 已 
在 智能 信息 处 理 、 医 学 、 工 程 、 地 质 灾害 测报 等 多 个 
领域 取得 了 良好 的 效果 (Roy et al., 2010; Elsafi， 
2014; Gosukonda et al.，2015 ) 。 小 波 分 析 (wavelet 
analysis) 具有 良好 的 时 频 局 部 特性 和 变焦 特性 , 且 
时 频 窗 可 根据 需要 及 时 调整 ,目前 已 成 功用 于 信和 号 
与 图 像 压缩 工程 技术 、 信 和 号 分 析 等 方面 (Fairley et 
al., 2014; Kosek et al., 2014)。Zhang 和 Benveniste 
(1992) 提 出 将 小 波 变换 理论 与 人 工 神经 网 络 相 结 
合 形 成 一 种 新 的 神经 网 络 , 把 小 波 函 数 作为 传递 函 
数 建立 小 波 神 经 网 络 (wavelet neural network， 
WNN) (Zhang and Benveniste, 1992; Bakshi and 
Stephanopoulos, 1993; Pati and 
1993 ) ,从 而 在 神经 网 络 研究 领域 发 展 了 一 个 新 的 
研究 方向 ,已 在 信号 处 理 、 图 像 压 缩 、 模 式 识 别 和 系 
统 辨 识 等 领域 开始 应 用 (Xu and Liu，2013; 
Guillermo et al., 2015 ) 。 

目前 ,人 工 神经 网 络 在 害虫 预测 预报 方面 已 有 
一 定 的 应 用 ( 陈 绘画 等 , 2003; 罗 盛 健 , 2006; 王国 
昌 , 2011; 陈 恩 会 等 , 2012) ,但 运用 小 波 神经 网 络 
对 预测 病虫害 发 生 情 况 在 国内 外 尚 处 于 起 步 阶 段 。 
本 研究 通过 建立 小 波 神经 网 络 模型 对 麦 蚜 发 生 程度 
进行 预测 ,并 与 BP 神经 网 络 模型 进行 对 比 ,对 两 者 
在 病虫害 测报 上 的 优 、 缺 点 进行 对 比分 析 , 以 期 为 更 
实用 的 预测 预报 模型 提供 依据 。 
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1 数据 来 源 


1.1 麦 蚜 数 据 

自 20 世纪 90 年代 中 后 期 以 来 ,从 发 生 面积 、 实 
际 损 失 来 看 ,小 麦 师 虫 已 上 升 为 我 国 农作物 重大 病 虫 
害 中 的 第 3 位 ( 曹 雅 忠 等 , 2006) ,在 山西 普 南 麦 蚜 也 
是 危害 最 严重 的 害虫 之 一 。 因 此 ,本 研究 将 麦 蚜 作 为 











研究 对 象 。 麦 蚜 原 始 数 据 来 自 山西 省 植保 植 检 总 站 。 

采样 地 点 :山西 省 运城 市 商城 县 古 魏 镇 北 关 村 ， 
34°36'N ~48°30'N,110°36'E ~ 42°30 玉 ,年 平均 气 
温 12.77% ,无 箱 期 250 d 左右 ,年 降水 量 513 mm。 
全 镇 耕地 面积 近 5 000 hm? ,地 势 平 坦 ,土壤 肥沃 ,为 
山西 小 麦 的 主 产 区 。 

历史 数据 的 调查 方法 : 单 对 角 线 5 点 取样 ,每 点 
国定 50 株 , 当 百 株 蚜 量 超过 500 头 时 ,每 点 可 减少 
至 20 株 , 定 点 、 定 时 及 固定 方法 调查 有 株 蚜 数 、 麦 蚜 
种 类 和 数量 。 

数据 范围 :本 研究 所 用 数据 自 1980 年 2 月 起 ， 
至 2014 年 6 月 止 , 于 每 年 2 月 底 到 6 月 初 每 隔 4 d 
采集 一 次 的 数据 。 

根据 农业 部 发 布 的 《小麦 蚜 虫 测报 调查 规范 》 
对 麦 蚜 原 始 数 据 统计 整理 ,将 研究 区 域 麦 蚜 发 生 程 
度 分 为 5 级 ,以 小 麦 蚜虫 发 生 盛 期 平均 百 株 小 麦 的 
蚜 量 (y) 来 确定 ,y<500 为 1 级 , 即 极 轻 度 发 生 级 ; 
500 <y 和 1 500 为 2 级 , 即 轻 度 发 生 级 ;1 500 < y= 
2 500 为 3 级 , 即 中 度 发 生 级 ;2 500 <y<3 500 为 4 
级 , 即 偏重 度 发 生 级 ;y >3 500 为 5 级 , 即 重 度 发 
生 级 。 

1.2 气象 数据 

气象 数据 来 自 山西 省 气象 局 。 以 1980 - 2014 
年 ,每 年 2 月 1 日 到 5 月 30 日 的 气象 数据 为 基础 , 统 
计 得 到 月 平均 温度 、 月 平均 最 高 温度 、 月 平均 最 低温 
度 月 平均 湿度 .月 平均 降水 量 月 平均 日 照 时 数 及 月 
平均 风速 等 单一 气象 指标 ;计算 复合 气象 指标 ,包括 
每 月 的 温 两 系数 (降雨 量 /平均 气温 )、 晴 十 系数 ( 降 
雨量 /日照 时 数 ) 和 温 湿 系数 (平均 湿度 /平均 气温 ) ; 
将 单一 气象 指标 与 复合 气象 指标 共同 建 表 ( 表 1) 。 

本 研究 将 1980 - 2009 年 麦 蚜 发 生 程度 作为 训 
练 集 ,2010 - 2014 年 麦 蚜 发 生 程度 作为 测试 集 , 具 
体 预测 流程 如 图 2 所 示 。 在 MATLAB 软件 平台 编 
写 神经 网 络 程序 ,其 建 模 程 序 分 为 数据 归 一 化 、 网 络 
训练 .网 络 预 测 、 误 差分 析 和 结果 作 图 等 过 程 ,采用 
紧 致 型 小 波 神经 网 络 3 层 结构 , 即 由 一 个 输入 层 一 
个 隐 含 层 和 一 个 输出 层 组 成 ,运用 LM 算法 进行 网 
络 训练 , 参 比 模型 为 BP 神经 网 络 。 





























































































































2 ”小波 神经 网 络 建 模 流程 


2.1 小 波 神经 网 络 简介 及 预测 流程 
小 波 分 析 于 1984 年 由 法 国 工程 师 英 莱 提 出 , 它 
能 有 效 地 在 时 间 和 空间 的 变换 中 提取 信息 ,在 伸缩 、 
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表 1 气象 因子 对 照 表 
Table 1 Comparison of meteorological factors 
平均 ”平均 最 高 气温 平均 最 低 气 温 平均 相对 湿度 平均 降水 量 平均 日 照 时 数 平均 风速 ” 温 雨 系数 ”晴雨 系数 ,yy 
2 , 温 湿 系数 
气温 Average Average Average Average Average Average Tempe- Sunshine 
月 份 有 . Temperature 
Average highest lowest relative Precipi- sunshine wind rature and rain a 
Month a Ne . .. and humidity 
tempera- temperature temperature humidity tation duration speed and rain coeffici- 放生 on 
ture (%C) (C) (C) (%) (mm) (h) (m/s) coefficient ent ey 
2 Xl 22 23 Xa Xs XG 37 Ng Xo X10 
3 X11 X12 EE: X14 2 X16 i X18 X19 X20 
4 X21 X22 X23 X24 X25 X26 X27 X28 X29 Xa0 
5 al 9g X33 X34 3535 X36 7 X38 X39 X40 


平移 等 运算 过 程 中 对 信号 进行 多 尺度 细 化 分 析 , 即 
在 时 间 窗 和 频率 窗 都 可 变化 的 固定 窗口 内 ,低频 部 
分 具有 和 较 罕 的 频率 窗 , 但 时 间 窗 较 宽 , 而 高 频 部 分 
反之 。 

基于 小 波 分 析 的 神经 网 络 有 两 种 结构 :一 种 是 
松散 型 小 波 神经 网 络 ( 图 1: A) , 即 先 用 小 波 分 析 对 
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图 1 





Fig. 1 


2.2 主 成 分 分 析 法 整合 气象 因子 

气象 因子 的 筛选 决定 着 病虫害 预测 模型 的 准确 
性 ,筛选 到 的 自 变 量 要 与 麦 蚜 的 发 生 具 有 较 高 的 相 
关 性 。 主 成 分 分 析 通 过 对 原始 多 维 输入 变量 形成 的 
数据 和 矩阵 求 取 相 关 和 矩阵 ,根据 相关 和 矩阵 的 特征 值 , 获 
得 累计 方差 贡献 率 ,再 根据 相关 和 矩阵 的 特征 向 量 , 确 
定 主 成 分 。 本 研究 中 ,作为 神经 网 络 的 输入 变量 的 
气象 因子 每 年 有 40 个 之 多 , 维 数 过 高 , 且 因 子 之 间 
具有 一 定 的 关联 性 , 若 将 它们 同时 作为 输入 变量 会 
造成 信息 的 过 分 使 用 ,从 而 增加 神经 网 络 模型 的 复 
杂 程 度 , 延 长 训练 过 程 ,降低 学 习 速 率 ,弱化 泛 化 性 
能 。 因 此 ,在 SPSS 软件 平台 运用 主 成 分 分 析 法 对 存 























信号 进行 预 处 理 , 然 后 通过 常规 神经 网 络 输出 ;为 一 
种 是 紧 致 型 小 波 神经 网 络 ( 图 1: B) ,是 目前 广泛 采 
用 的 一 种 结构 形式 , 即 用 小 波 函 数 代替 常规 神经 网 
络 中 的 隐 含 层 函 数 , 用 小 波 函 数 的 尺度 代 蔡 输入 层 
到 隐 含 层 的 权 值 , 用 小 波 函 数 的 平移 参数 代替 隐 层 
国 值 ( 传 苓 璇 等 , 2010) 。 











松散 型 小 波 ( A) 和 紧 致 型 (B) 人 小波 神经 网 络 


Loose-type (A) and compact-type (B) wavelet neural network 


一 化 处 理 到 [0,1] 范围 内 ( 罗 长 寿 等 , 2004; 陈 明 ， 
2013 ) 。 归 一 化 公式 为 : 


式 中 ,x; 表示 数据 原始 值 ,x 表示 归 一 化 后 的 数据 ， 
xs 和 win 分 别 表示 每 一 类 成 分 因子 的 最 大 值 和 最 
小 值 。 
2.4 ” 试 凑 法 筛选 神经 网 络 传递 函数 和 隐 含 层 节 
点 数 

神经 网 络 的 参数 很 多 ,有 输入 层 节 点 数 、 隐 含 层 
节点 数 、 输 出 层 节 点 数 \ 传 递 函数 训练 函数 及 学 习 
率 等 (Liu et al., 2006) 。 在 网 络 结构 .权重 及 阔 值 相 


多 





在 共 线 性 的 气象 因子 重新 整合 ,形成 新 的 成 分 因子 
后 输入 模型 ( 刘 了 凡 等 ,1997 ) 。 

2.3 ” 归 一 化 法 处 理 数据 

由 于 不 同 成 分 因子 量 纲 不 同 ,神经 网 络 对 [0,1] 
间 的 数据 最 敏感 ,因此 ,在 建 模 之 前 ,将 成 分 因子 归 














同 的 情况 下 ,选择 适当 的 隐 含 层 节 点 数 及 传递 函数 ， 
直接 影响 着 网 络 的 学 习 能 力 和 泛 化 能 力 。 本 研究 运 
用 试 凑 法 对 神经 网 络 的 隐 含 层 节点 数 和 传递 函数 进 
行 参 数 寻 优 ,以 期 得 到 最 好 的 初始 值 组 合 ,其 余 参数 
根据 经 验 取 值 。 
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i 表 2 主 成 分 分 析 的 累积 百分率 
Input original data Table 2 Cumulative percentage of principal 
component analysis (PCA) 
> 页 献 率 ( % 累积 % 
主 成 分 分 析 法 整合 形成 新 成 分 因子 ee he 
Variance Cumulative variance 
Create new factor by PCA 
1 22.632 22.632 
2 14.772 37.404 
归 一 化 法 处 理 数据 3 12. 558 49.962 
Manipulate data through normalization method 4 9 552 59 514 
S$ 9.020 68.534 
采用 试 凑 法 寻找 最 优 参数 6 7.683 76.217 
Select the optimal parameters by 

trial and error method 4 人 
8 3.932 84.756 
9 2.869 87.625 









是 否 找到 最 优 参数 
Whether successfully select 
the optimal parameters 





建立 小 波 神经 网 络 
Build up wavelet neural network (WNN) 


对 测试 集 进行 预测 


Forecast the test set 


输出 预测 结 
Output result of the forecast 


图 2 小 波 神 经 网 络 预测 流程 
Fig. 2 Prediction flow of wavelet neural network (WNN) 








2.5 模型 验证 方法 
本 研究 选择 平均 绝对 百 分 误 差 (MAPE ) 和 均 方 
误差 (MSE) 两 种 评价 指标 对 模型 性 能 进行 评价 (向 

昌盛 和 周子 英 , 2010; 李 启 权 等 , 2013 ) 。 
MAPE = Ley A 


ni=1 yi 











MSE = Dy (y; 7) 。 
公式 中 , y% 为 发 生 程度 的 实际 值 , 六 为 模型 的 预测 
值 ,n 为 训练 样本 数 。MAPE 值 越 小 ,说 明 模 型 预测 
精度 越 高 ;MSE 值 越 小 ,说 明 模 型 稳定 性 越 好 。 











3 结果 


3.1 主 成 分 分 析 结 果 

将 表 1 气象 因子 对 照 表 中 的 40 个 自 变量 进行 
主 成 分 分 析 , 得 到 主 成 分 分 析 的 累计 百分率 ( 表 2)， 
前 9 个 成 分 的 累计 贡献 率 达 到 87. 625% ,这 9 个 新 











成 分 不 仅 基本 包含 全 部 影响 因素 的 信息 ,而且 彼此 
之 间 互 不 相关 。 其 余 31 个 主 成 分 的 贡献 率 只 有 
12.375% ,可 作为 噪声 项 不 进入 建 模 过 程 。 表 3 为 
按照 成 份 矩阵 重新 形成 的 成 分 因子 z1 ~z 与 麦 蚜 发 
生 程 度 组 成 的 新 的 样本 集 。 
由 表 3 成 份 矩 阵 表 可 知 每 个 新 成 分 因子 与 原 有 
因子 之 间 的 相关 系数 及 新 因子 形成 的 表达 式 。 如 
zi! =X1 * 0.008 — x, *0.069 +xs *0.055... Nag * 
0.400 +x3o 0.371 +xao * 0.706, 以 此 类 推 ,分 别 得 
到 z ~z 的 表达 式 。 同 时 ,相关 系数 绝对 值 最 大 的 
前 3 位 分 别 是 ,zi 主 成 分 中 为 x1,(3 月 平均 最 高 温 
度 ) xx(3 月 温 湿 系 数 ) x3(5 月 平均 最 高 温度 ) ;2 
主 成 分 中 为 Xao(5 月 晴雨 系数 ) \x35(5 月 平均 降水 
量 ) ws(5 月 温 雨 系数 ) ;zs 主 成 分 中 为 (4 月 睛 
雨 系数 ) .xs (4 月 温 雨 系数 ) .xo (2 月 晴雨 系数 ) ;有 
主 成 分 中 为 x1(2 月 平均 气温 ) x,(2 月 平均 最 高 气 
温 ) xs(2 月 平均 最 低 气温 ) ;zs 主 成 分 中 为 x。(2 月 
晴雨 系数 ) x (4 月 平均 最 低 气温 ) x;(2 月 平均 降 
水 ) ;ze 主 成 分 中 为 x (4 月 平均 气温 ) .xss (5 月 平 
均 最 低 气温 ) ,xz(3 月 平均 风速 ) ;zj 主 成 分 中 为 xy 
(5 月 平均 风速 ) xio(2 月 温 湿 系数 ) x (2 月 平均 
风速 ) ;zs 主 成 分 中 为 za (5 月 平均 最 低 气温 ) was 
(3 月 平均 日 照 时 数 ) x (3 月 平均 风速 ) ;zo 主 成 分 
中 为 xio(2 月 温 湿 系 数 ) x (2 月 平均 风速 ) x3 (5 
月 平均 相对 湿度 ) 。 

从 时 间 上 来 看 ,2 和 5 月 气象 因素 对 麦 蚜 发 生 
的 影响 较 3 和 4 月 大 ,2 月 天 气 逐 渐 转 暖 ,越冬 麦 蚜 
开始 活动 ,繁殖 机 会 增加 , 可 为 穗 蚜 积累 更 多 的 虫 
源 ,5 月 达到 麦 蚜 生 长 发 育 的 最 适 条 件 , 且 小 麦 处 于 
抽穗 ,扬花 和 灌浆 的 关键 时 期 , 麦 蚜 密 度 直 线 上 升 ; 
从 因素 上 看 ,与 温度 湿度 相关 的 气象 因子 对 麦 蚜 发 
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表 3 成 分 矩阵 表 
Table 3 Component matrix 
Z1 22 23 Z4 Zs Z6 27 28 Zo 
Xl 0.008 0.367 —0.009 0.844 0.027 —0.140 —0.246 0.032 0.039 
2 一 0.069 0.262 —0.233 0.779 0.194 一 0. 190 一 0.287 0. 060 0.035 
NX3 0.055 0.435 0.256 0.716 —0.253 —0.072 —0.126 0.008 0.014 
4 —0.025 0.025 0.627 0.061 —0.502 0.168 0. 174 0. 136 —0.258 
Xs —0.038 0.001 0.621 0.207 —0.547 0.171 —0.137 —0.134 —0.245 
X6 一 0.039 一 0.252 一 0.484 0. 149 0.537 —0.269 —0.240 0.088 —0.080 
NX7 0.273 —0.036 -0.110 —0.264 一 0.425 —0.181 0.444 —0.179 0.339 
X8 0.339 —0.106 0.549 0.012 —0.512 0.180 —0.319 —0.131 0.201 
Xo 0.014 0.092 0.661 0. 101 —0.582 0.206 —0.057 —0.172 —0.100 
X10 0.387 —0.165 0.185 0.167 —0.226 0.072 —0.514 0.144 0.477 
X11 —0.746 0.231 0.303 0.312 0. 123 -0.310 0. 123 0.023 —0.055 
X12 —0.804 0.234 0.277 0.204 0.086 —0.285 0.168 —0.046 一 0.037 
X13 —0.564 0.256 0.223 0.481 0.201 —0.245 0.065 0.140 —0.043 
X14 0.694 —0.236 —0.183 0.345 0.031 0.377 —0.003 0.065 一 0.204 
X15 0.592 —0.244 —0.252 0.560 0.019 0.154 0.211 —0.157 —0.004 
X16 —0.736 0.060 0.345 —0.137 0.072 一 0.034 0. 099 一 0.437 0.014 
X17 0.347 —0.238 —0.175 —0.151 —0.124 —0.473 0.266 0.414 0.164 
X18 0.691 —0.285 —0.337 0.393 —0.041 0.223 0.131 —0.134 -0.010 
X19 0.632 —0.245 —0.342 0.517 0.003 0.174 0.169 —0.088 —0.031 
X20 0.777 —0.263 —0.323 —0.062 —0.095 0.382 —0.090 0.069 —0.053 
Nol —0.396 0.566 —0.251 0.031 0.242 0.555 0. 104 一 0.204 0. 069 
X27 —0.505 0.539 —0.350 0.117 0.060 0.465 0.047 —0.113 0.069 
X23 —0.030 0.485 0.071 一 0.029 0.560 0.448 0.179 —0.304 0.078 
X24 0.397 —0.368 0.429 0.171 0.486 0.024 0.248 —0.053 —0.267 
N25 0.268 —0.045 0.658 —0.002 0.525 0.331 0.045 0.143 0.231 
X26 -0.233 0.174 —0.509 0.069 —0.257 0.376 0.049 —0.167 0.142 
X27 0.555 0.171 —0.004 0.151 —0.165 —0.435 0.371 —0.172 0.106 
X28 0.292 一 0. 104 0.687 —0.028 0.481 0.287 0.029 0.173 0.217 
X29 0.272 —0.081 0.701 —0.043 0.504 0.293 0.036 0.165 0.179 
X30 0.476 —0.550 0.426 0.132 0.235 —0.252 0.155 0.014 一 0.242 
X31 —0.738 —0.388 —0.122 0.082 —0.060 0.332 —0.003 0.302 —0.052 
X32 —0.777 —0.448 —0.081 0.179 —0.035 0.231 —0.048 0.179 0.004 
X33 —0.389 0.139 —0.226 —0.233 —0.200 0.514 0.050 0.512 —0.270 
X34 0.618 0.502 —0.003 -0.310 0. 161 0. 022 —0.231 -0.130 -0.313 
X35 0.366 0.845 0.057 —0.009 —0.011 —0.033 0. 139 0.237 0.037 
X36 —0.451 —0.671 0.107 0.203 0.148 0.059 0.125 一 0.227 0. 095 
X37 —0.043 —0.010 0.008 0.406 -0.315 0.266 0.598 0.187 0.056 
X38 0.400 0.845 0.066 —0.001 0.010 —0.045 0. 138 0.201 0.036 
X39 0.371 0.847 0.053 —0.020 —0.034 —0.038 0.121 0.262 0.022 
X40 0.706 0.511 0.044 —0.232 0.148 —0.105 —0.161 -0.206 一 0.203 





生 的 影响 较 大 ,如 气温 、 湿 度 、 降 水 及 温 湿 系数 、 温 雨 
系数 等 。 同 时 ,由 于 麦 蚜 是 一 种 迁 飞 性 害虫, 风速 对 
麦 蚜 的 生长 发 育 有 一 定 的 影响 ;日 照 时 数 在 一 定 范 
转 内 ,可 以 促进 雌 蚜 的 繁殖 。 

表 4 为 按照 成 份 矩 阵 重新 形成 的 成 分 因子 za ~ 
z 与 麦 蚜 发 生 程 度 组 成 的 新 的 样本 集 。 
3.2 ”模型 参数 筛选 结果 
3.2.1 小 波 神经 网 络 参 数 筛 选 结 果 :小 波 神经 网 络 


























和 BP 神经 网 络 最 大 的 区 别 在 于 传递 函数 的 不 同 。 
小 波 神经 网 络 选择 Morlet 小 波 函 数 作为 传递 函数 
(图 3) 。 

它 是 高 斯 函数 下 的 单 频率 副 正弦 函数 ( 陈 哲 和 
汉 天 弄 ， 1999) ,C 是 重 构 时 的 归 一 化 常数 ,其 定义 
如 下 : 








w(t) = Cb 008 C3 
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表 4 ， 主 成 分 分 析 后 的 麦 蚜 发 生 程度 样本 集 
Table 4 Sample sets of the occurrence degree of wheat aphids after principal component analysis (PCA) 

年 份 Year 21 22 23 24 25 26 2 28 29 y 
1980 1.3570 —1.2524 —0.2601 1.5610 —0.0542 —2.5811 1.4322 —1.3436 —0.0718 4 
1981 0.6621 —1.0246 -1.0006 -0.7995 —0.1180 一 0.8396 0.5417 1.7544 1.4291 1 
1982 一 1.0792 —1.4786 0.2606 0.4980 —0.0664 一 1.1747 0.8357 —0.3759 1.2747 2 
1983 0.0974 —0.7690 —0.8672 0.5974 —1.1470 —0.1284 0.9272 一 0.2872 0.4933 3 
1984 0.7800 0.1987 一 0.8981 一 0.7069 0. 6337 -0.9400 -0.0990 —1.2745 0.9772 1 
1985 —0.0586 0.3285 0.0102 -2.3099 1.0828 一 1.0479 1.2173 —0.5856 一 0.8820 1 
1986 0.9475 1.3445 -1.1295 -1.7087 一 0.9507 0.5152 -0.3113 0.5230 0.4254 1 
1987 —0.1762 —0.6906 -0.3427 -1.1990 0.7726 -0.8592 -1.3210 0.0658 0.3524 1 
1988 0.9462 —0.0120 -0.5045 -0.2507 —0.0172 -0.1432 -0.3738 0.2231 一 0.7972 1 
1989 1.3252 —0.5198 -0.7967 -0.8644 —1.3359 0.9365 —0.8133 0.0922 —0.2653 1 
1990 0.8687 一 0.3599 1.6067 -0.9814 0.2401 0.7701 —0.6293 —0.8951 0.4741 1 
1991 0.8498 —0.4458 1.9917 一 0.0673 —2.2935 0.1499 —0.1307 —0.3293 —0.8857 1 
1992 1.5651 —0.7756 —0.5803 0.9643 —0.4383 一 0.2019 0.9267 —0.7548 —0.5369 1 
1993 1.0617 —0.6436 —1.7261 0.6761 0.7122 1.4238 —0.5818 0.2843 —0.7813 2 
1994 1.0128 0.7812 0.9208 —0.0643 0.3033 —0.6059 —0.6436 —0.8436 —0.2505 1 
1995 —0.0324 —1.6881 1.1645 0.2445 0.7486 2. 1068 0.8168 0.5689 1.2446 4 
1996 0.1564 —1.0601 0.5269 -1.2372 1.0748 一 0.2299 -1.6217 1.0939 —0.1877 3 
1997 —0.3083 一 0.9952 0.3060 0. 8960 —0.6143 0.5740 -0.5268 —0.3657 —2.1794 2 
1998 1.5474 1.9761 0.6114 0. 8060 1.8372 0.9961 1.4582 0.5682 0. 1210 1 
1999 0.3063 0.7522 一 0.4717 1.3243 2.2627 0.5568 -0.5800 0.0676 0.5680 1 
2000 —1.9112 一 0.4720 —0.0489 -0.3214 0.5293 0.1320 -0.7260 —0.0051 0. 8855 2 
2001 —1.6049 —0.6126 0.5648 0.2869 —0.1457 —0.9648 -1.1585 —0.6057 —1.0826 9 
2002 0.0280 1.9041 一 0.0474 0.8833 0.5375 —1.8011 —0.7856 —0.3761 一 0.3964 2 
2003 0.5130 一 0.0947 1.5706 0.6365 0.0747 0.7036 -0.3984 1.2309 —1.5464 S| 
2004 —1.2471 0.8369 —0.6198 1.4407 0.7984 0.2652 -1.3409 —1.0075 0. 8821 2 
2005 —1.6803 0.6610 -1.3444 -1.3161 一 0.3274 1.0834 2. 8002 一 0.4593 —1.3983 2 
2006 一 0.9992 0.8358 0.7623 -0.2027 0.0130 0.3414 -0.6326 —0.5618 一 0.3771 3 
2007 -1.1513 —0.2936 —1.4223 2.3054 —1.1161 0.9009 —0.3000 0.8410 0. 1208 2 
2008 —1.6123 —0.5387 -0.1480 -0.6230 1.3947 -0.1130 0.7734 0.1198 —2.2300 2 
2009 0.0478 2.3007 0.5267 0.6595 —1.0633 0.0271 0.5486 0.5780 —0.2679 2 
2010 0.0268 —0.0323 0.4263 —0.3159 —0.5197 一 0.9953 —0.6831 2.2730 0.2936 2 
2011 —0.5682 0.9726 -0.0502 -0.7387 一 1.7809 0.9760 -0.8392 —2.5767 1.3851 2 
2012 0.0142 —0.2483 —1.3850 -0.3462 一 0.0284 1.0391 0. 0429 0. 1383 0.7210 3 
2013 一 1.0943 1.3142 0.1885 0.1579 —1.3665 一 1.4117 0.5014 2.4194 0.6798 3 
2014 —0.5898 —0.1990 2.4053 0.1160 0.3674 0.5397 1.6743 —0.1939 1.8091 2 

Morlet 时 域 Morlet 频 域 
Morlet time domain Morlet frequency domain 
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图 3 Morlet 函数 


Fig. 3 Morlet function 
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MATLAB 的 传递 函数 中 没有 Morlet 小 波 函 数 ， 
需要 创建 自 定义 的 传递 函数 ,编写 的 2 个 子 程序 分 
别 是 : 

(1)function y =mymorlet(1); 

y=exp( -(t.2)/2) * cos(1.75 *1)。 

(2)function y = d_mymorlet(1); 

y= —1.75 *sin(1.75 *1). *exp( -(t.2)/ 
2) -tx*cos(1.75*1). *exp( -(t.2)/2)。 
增加 隐 含 层 节 点 数 能 有 效 提高 神经 网 络 的 训练 
精度 ,但 隐 含 层 节 点 数 过 多 又 会 降低 训练 速率 。 选 
取 隐 含 层 神经 元 节点 数 时 参考 以 下 公式 : 

(1) 1 = | rz| ( 李 伟 超 等 , 2006) ; 

(2) 1 < | Vmr+a|+a( 陈 志高 , 2009 ) ; 

(3)! < n--1( 史 峰 等 , 2010); 

(4)71 = 2n+1( 王 德 明 等 , 2012 ) ; 

(5)1 = log2n ( 王 尚 昆 , 2013 ) 。 

式 中 ,n =9( 输 入 层 节 点 数 ),m =1( 输 出 层 节 点 
数 ) ,a 在 [1,10] 之 间 , 计 算 不 同 隐 含 层 节 点 数 对 应 
的 预测 误差 。 

小 波 神经 网 络 的 隐 含 层 节 点 数 在 公式 计算 的 基 
础 上 ,根据 经 验 进行 调整 ,对 隐 含 层 节 点 数 在 10 ~ 
20 之 间 进 行 测试 ,结果 显示 , 当 隐 含 层 节点 数 为 15 
时 ,小 波 神 经 网 络 训练 的 MSE 和 MAPE 值 均 达到 最 
小 ( 表 5)。 

经 过 第 选 和 经 验 取 值 ,建立 了 3 层 的 小 波 神经 





























y=logdig(s) 


y=tansig(s) 


网 络 结 构 9-15-1 ,学 习 概 率 1r1 =0. 01 ,12 =0.05 ;和 迭 
代 次 数 150 ,误差 累计 为 0。 

表 5 ”小 波 神经 网 络 中 不 同 隐 含 层 节点 数 对 应 的 预测 误差 
TableS Forecast errors corresponding to different hidden 


layer nodes in wavelet neural network ( WNN) 


隐 含 层 节 点 数 








Node number of hidden layer neurons Wh WALE (CW) 
10 0. 5200 22. 8955 
11 1.2321 25.7885 
12 0.9904 34. 8834 
13 0. 5786 24. 9330 
14 1.0774 26. 2084 
15 0.0578 10. 1729 
16 0. 8960 27.9961 
17 0.9607 33.7603 
18 0. 9455 34. 1062 
19 1.0125 29.7315 
20 0.9302 34.3630 


3.2.2 BP 神经 网 络 参 数 筛选 结果 : BP 神经 网 络 
一 般 选 用 Sigmoid 作为 传递 函数 , 它 主要 有 3 种 类 型 
(图 4) ,分 别 为 logsig 函数 、tansig 函数 及 purelin 

Sigmoid 六 数 是 光滑 、 可 微 的 , 比 线性 函数 更 精 
确 ,容错 性 更 好 ,其 定义 如 下 : 


1 
S = 
(4) 1 +e” 








t 站 


y 





y=purelin(s) 


图 4 Sigmoid 函数 


Fig. 4 Sigmoid function 


将 3 种 函数 循环 组 合 , 得 到 不 同 传递 函数 对 应 
的 预测 误差 ( 表 6) ,虽然 logsig-tansig、purelin-logsig 
这 两 个 组 合 的 MSE 值 较 tansig-purelin 小 ,但 MAPE 
值 却 大 于 后 者 ,因此 选用 tansig-purelin 作为 传递 函 
数组 合 。 

对 BP 神 经 网 络 隐 含 层 节 点 数 进行 筛选 , 当 节 
点 数 为 10 时 ,BP 神经 网 络 训练 的 MSE 和 MAPE 值 
均 达 到 最 小 ( 表 7)。 

经 过 筛选 ,确定 BP 神经 网 络 结构 为 9-10-1, 隐 








含 层 传递 函数 为 tansig, 输 出 层 传递 函数 purelin , 训 
练 图 数 为 rainlm ,网 络 的 权 值 学 习 涌 数 为 learngdm， 
性 能 函数 为 MSE ,动量 因子 为 0.9; 训 练 次 数 为 30; 
学 习 率 lr =0.1; 期 望 目标 误差 最 小 值 1 x10“。 
3.3 1980 -2009 年 麦 蚜 发 生 程 度 拟 合 结果 

运用 小 波 神经 网 络 的 训练 过 程 中 ,有 21 年 拟 合 
结果 的 误差 实际 值 在 0 ~ 1 之 间 , 其 中 11 年 拟 合 精 
度 超过 80% (图 5: A)。 运 用 BP 神经 网 络 对 1980 - 
2009 年 麦 蚜 发 生 程度 进行 训练 (图 5: B) ,有 25 年 
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表 6 不 同 传递 函数 对 应 的 预测 误差 表 7 BP 神经 网 络 中 不 同 隐 含 层 节 点 数 对 应 的 预测 误差 
Table 6 Forecast errors corresponding to Table 7 Forecast errors corresponding to different hidden 
different transfer function layer nodes in BP neural network (BPNN) 
总 含 层 传递 函数 ”输出 层 传递 函 妆 颖 含 层 节点 类 MAPE 
序号 隐 含 县 传递 函数 输出 县 传递 函数 这 1 区 含 层 闻 数 i 
N Transfer function Transfer function MSE (% ) Node number of hidden layer neurons (%) 
of hidden layer of output layer . 1 0.9971 36.7857 
1 logsig tansig 0.5447 26.2215 2 2.4013 50.0235 
2 logsig purelin 0.7599 27.1611 3 1.3537 35.0445 
3 logsig logsig 0.7629 33.8816 4 0 7437 23 9774 
4 tansig tansig 1.2431 47.3269 
5 0.6407 33.8786 
5 tansig logsig 1.4327 47.5664 
6 tansig purelin 0.6192 16.9325 ® 1 23:780 
7 purelin logsig 0.6000 ”30.0000 1.5801 32. 6757 
8 purelin tansig 1.0464 40.4621 8 2.4714 43.6849 
9 purelin purelin 0.7981 28.8322 9 0.7027 28.5512 
10 0.6192 16.9325 
的 发 生 程度 拟 合 结果 的 误差 实际 值 在 0 ~1 之 间 ( 即 11 0.7743 28.4877 
误差 在 1 级 以 内 ) ,其 中 13 年 拟 合 精度 超过 80% 。 说 人 人 nk 
He A ,不 、 < 13 1.5889 48.3045 
明 两 个 模型 都 能 较 好 描述 麦 蚜 发 生 程度 的 变化 规律 ， ~ es 有 


BP 神经 网 络 的 拟 合 能 力 较 小 波 神经 网 络 好 。 
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图 5 小 波 神经 网 络 (A) 和 BP 神经 网 络 (B) 训练 拟 合 结果 图 
Fig. 5 Simulation effect of training of wavelet neural network (WNN) (A) and BP neural network (BPNN) (B) 








3.4 2010 -2014 年 麦 蚜 发 生 程度 预测 结果 (1) 两 种 模型 的 平均 预测 精度 都 超过 80% ,小 
比较 小 波 神经 网 络 和 BP 神经 网 络 的 预测 精度 ”小 神经 网 络 平均 预测 精度 为 89. 83% ,BP 神经 网 络 
和 性 能 可 知 ( 表 8 和 9) : 的 平均 预测 精度 为 83.07% 。BP 神经 网 络 预测 的 
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2010, 2011 及 2014 年 的 麦 蚜 发 生 程度 较为 准确 , 均 
超过 90% ,但 2013 年 的 预测 值 与 实际 值 相差 较 大 ， 
拟 合 精度 过 低 , 而 导致 平均 预测 精度 小 于 小 波 神 经 
网 络 的 平均 值 ,这 是 由 于 BP 神经 网 络 易 于 出 现 过 








致 的 预测 偏差 。 

(2) 比较 两 种 模型 的 均 方 误差 (MSE) ,小 波 神 
经 网 络 的 MSE 值 小 于 BP 神经 网 络 , 其 预测 的 均 好 
性 高 于 BP 神经 网 络 , 表 明 其 稳定 性 优 于 BP 神经 











度 训 练 和 训练 不 足 、 易 于 陷 人 局 部 最 小 等 缺陷 而 导 


网 络 。 


表 8 小 波 神经 网 络 和 BP 神经 网 络 预测 拟 合 结果 
Table 8 Fitted result of wavelet neural network (WNN) and BP neural network (BPNN) forecast 











小 波 神经 网 络 BP 神经 网 络 

年 从 wavelet neural network ( WNN) BP neural network ( BPNN) 
二 全 有 了 本 A , a i BE 
Year 实际 值 。 预测 值 。 误差 绝对 值 。 拟 合 精度 (%) | 汪汪。 预测 值 。 。 误 关 绝对 值 。 拟 合 精度 (9 ) 

Actual Predicted Error absolute Fitting 人 Predicted Error absolute Fitting 

Actual value 

value value value accuracy value value accuracy 
2010 2 2:32 0.32 84.18 2 2.05 0.05 97527 
2011 2 1.98 0.02 98.95 2 2::17 0.17 91.60 
2012 3 2:55 0.45 84. 84 3 2.50 0.50 83. 46 
2013 3 2.59 0.41 86.20 3 1:32 1.68 44.05 
2014 2 2.10 0.10 94.99 2 2502: 0.02 98.96 


表 9 BP 神经 网 络 和 小 波 神 经 网 络 性 能 比较 
Table 9 Performance comparison of BP neural network 
(BPNN) and wavelet neural network (WNN) 








立 了 9-15-1 的 小 波 神经 网 络 和 9-10-1 的 BP 神经 网 
络 ,从 结构 上 说 ,小 波 神 经 网 络 的 神经 元 为 15 个 ， 
BP 神经 网 络 的 神经 元 为 10 个 ,前 者 的 结构 更 为 复 
杂 , 说 明 对 数据 的 分 析 更 为 精确 ; 三 是 筛选 传递 函 








模型 Model MAPE (% ) MSE 
PCA-BPNN 16.93 0.6192 
PCA-WNN 10.17 0.0578 


4 ”结论 与 讨论 








农作物 病虫害 的 发 生 具 有 复杂 的 非 线性 特点 ， 
本 研究 基于 人 工 神经 网 络 和 人 小波 分 析 理 论 相 结合 ， 
提出 了 融合 两 者 而 形成 的 一 种 非 线 性 预测 模型 一 一 
小 波 神 经 网 络 模型 。 本 研究 充分 利用 山西 省 麦 蚜 发 
生 监 测 数据 建立 了 小 波 神经 网 络 模型 与 BP 神经 网 
络 模型 ,并 进行 了 预测 精度 和 稳定 性 对 比 。 结 果 显 
示 , 前 者 在 预测 准确 率 和 模型 稳定 性 上 均 好 于 后 者 ， 
其 同时 具有 小波 分 析 时 频 同步 分 析 的 特点 并 兼 具 神 
经 网 络 良好 的 自学 习 、 自 组 织 及 非 线性 映射 能 力 ,该 
模型 可 作为 农作物 病虫害 预测 预报 的 一 种 新 方法 继 
续 探 索 尝 试 。 

本 研究 中 小 波 神 经 网 络 模 型 能 够 较 好 预测 麦 蚜 
发 生 程度 的 本 质 在 于 :一 是 在 建立 神经 网 络 前 ,对 输 
入 变量 运用 主 成 分 分 析 法 进行 预 处 理 ,不 仅 将 有 一 
定 关联 的 因子 之 间 的 相互 影响 降 到 最 小 , 而 且 使 维 
数 由 40 个 降低 到 9 个 ,避免 了 由 于 “高 维 数 ” 出 现 
的 “ 维 数 灾 ”; 二 是 筛选 隐 含 层 节 点 数 ,由 于 在 理论 
上 尚 无 准确 且 完 备 的 确定 隐 含 层 节 点 数 的 方法 ,在 
已 有 研究 基础 上 ,运用 试 凑 法 筛选 出 节点 数 ,分 别 建 










































































数 ,小 波 神经 网 络 选择 常用 的 小 波 函 数 Morlet 函数 
作为 传递 函数 ,BP 神经 网 络 选用 tansig 作为 隐 含 层 
传递 函数 ,purelin 作为 输出 层 传递 函数 。 前 者 使 用 
的 Morlet 小 波 函 数 比 后 者 使 用 的 Sigmoid 函数 对 非 
线性 特征 的 提取 能 力 更 强 , 对 干扰 的 鲁 棒 性 好 ,误差 
小 , 泛 化 能 力 好 。 

然而 ,将 小 波 神经 网 络 运用 于 病虫害 预测 预报 
还 存在 许多 或 待 解 决 的 问题 。 例 如 :如 何 更 准确 地 
确定 神经 网 络 隐 含 层 节 点 个 数 ;在 预测 不 同 的 指标 
时 选择 何 种 小 波 函 数 ;初始 权 值 .尺度 伸缩 因子 、 时 
间 平 移 因 子 等 参数 的 选择 , 以 使 病虫害 预测 预报 更 
加 精准 和 稳定 等 ,希望 在 以 后 的 研究 中 进一步 探索 。 
研究 中 没有 考虑 麦 蚜 越冬 基数 天敌 寄主 及 防治 措 
施 等 因素 的 影响 ,可 通过 进一步 实验 对 这 些 因素 做 
定性 定量 分 析 后 进入 建 模 过 程 。 
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